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RÉSUMÉ* - Une nouvelle méthode est proposée pour estimer par itération le paramètre 
de la courbe de croissance de Von Bertalanffy Elle est basée sur la résolution d une équation 
de la forme L w ) par approximations successives. Elle a été appliquée aux Trachurus 

trachurus (Carangidae) de Tunisie, 


ABSTRACT. - A new method for fish grovvth parameters estimation. 

A new method is proposed in order to estimate by itération the parameter of the 
Von Bertalanffys growth curve. The method is based on the resolution of an équation of the 
shape by successive approximations. Il has been applied to Trachurus trachurus 

(Carangidae) from Tunisia, 


Mots-clés. - Carangidae, Trachurus trachurus, MEB, Tunisie, Growth, Von Bertaîanfïy's pa¬ 
rameters. 


La courbe de Von BertalanfTy (1938; 1949) est le modèle le plus utilisé par 
les ichtyologistes pour représenter îa croissance linéaire des Poissons, notamment 
en dynamique des populations. De nombreuses méthodes ont été mises au point 
pour l'estimation des trois paramètres de cette courbe. Nous citons à titre 
d'exemple les suivantes qui sont d'ailleurs les plus utilisées: 

1- Les méthodes graphiques introduites par Beverton et Holt (1957). k 
et t c sont estimés graphiquement en traçant a vue la droite d'ajustement qui 
représente au mieux un ensemble de points dont les coordonnées sont déduites des 
données d'observation, â savoir les longueurs moyennes à des âges régulièrement 
espacés, 

2- La méthode de Tomlînson et Abramson (1961) qui utilise la méihode des 
moindres carrés pour une estimation conjointe des trois paramètres à partir des 
données individuelles, sans passer par les longueurs moyennes atteintes à des âges 
régulièrement espacés les uns des autres. 

3- Les méthodes développées par Gulland (1965) qui consistent â calculer 
les valeurs des trois paramètres à partir de l'équation des droites d'ajustement 
déterminées par la méihode des moindres carrés (droites de régression), 

4- La méthode de Rafail (1973) qui consiste à déterminer les paramètres k 
et L w par itérations successives â partir d'une première estimation de ces 
paramètres, 

La meilleure méthode est celle de Tomlînson et Abramson qui fait appel au 
principe du maximum de vraisemblance. Elle présente toutefois l'inconvénient de 
nécessiter des calculs relativement compliqués. Celle que nous proposons ici est 
plus simple et conduit à des résultats pratiquement équivalents, A l'inverse de la 


(Il Biologie animale, Faculté des Sciences de Tunis, 1060, Campus universitaire, Tunis, TU¬ 
NISIE. 

(2) Muséum national d Histoire naturelle de Paris, 43 rue Cuvier, 75231 Paris Cedex 05, 
FRANCE. 


Cybium !991 t 15(2): 103-109 , 



104 


méthode de Rafail qui utilise t \ussi un processus d'itération, notre méthode est fa¬ 
cilement applicable à tous les cas. 


Rappel des formules de base 

L'équation de la courbe de Von Berralanffy est de la forme : 

L,~ LJi l0] (f) 

En dérivant les deux membres de celte équation par rapport à t et en prenant les 
logarithmes népériens, on obtient facilement: 

log-^ = lQg L x + log A — fcf + ki 0 

On en déduit: 


, l'n + l ^"n *n - 1 + *n 

log-j-—7~ = «-2- + b 

‘rc ï-1 z 


( 2 ) 


k — — a est donc ta pente de la droite de régression du logarithme de 
1 accroissement de la longeur entre deux âges successifs par rapport à l age moyen 
correspondant. 

De Téqualion (t) on tire egalement: 

log(l (3) 


k est donc la pente de la droite de régression du logarithme népérien de 

( . L t \ 

f 1 — j —- I par rapport a L 

Enfin lorsqu'on connaît une valeur de k et donc de e k t on peut calculer la 
valeur correspondante de L^ par la formule suivante déduite de (1): 


L 


oc 




«* 11 ,.:- Z L, 

et pour un échantillon L x — -— - ——-* - 1 t — 

OU 


(B —!)(«*- 1) 


(4) 


I) — 1 71 — 1 

en posant S\ — Y. L t + 1 et S 2 — Y. L t 
1 1 

Méthode de Rafail (1973) 

Celle-ci consiste à calculer ur s première valeur de k par la relation (2), puis 
une première valeur de L x par la relation (4). On calcule ensuite la valeur corres¬ 
pondante de t 0 ainsi que V*/ 2 , somme des carrés des écarts entre les valeurs L t 
observées et les valeurs L t données par le modèle pour les memes âges. Cette 
somme permet d'évaluer globalement le degré d'ajustement du modèle aux 
données de l'observation. 

Rafail calcule une seconde valeur de k par la relation (3), puis, comme 
précédemment, une seconde valeur pour f ffl cl V/^ 2 , etc. S'il y a convergence 
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"£>d 2 se stabilise à une valeur minimale au bout de plusieurs itérations. Les esti¬ 
mations de L mî k et l 0 sont alors, elles aussi, stabilisées et le processus d'itération 
peut être arrêté. 

Rafail a traité deux exemples relatifs à Serranus alexandrinus et flippoglos- 
sus stenolepis où effectivement diminue à chaque itération, les calculs étant 
arrêtés après la seconde ou la troisième. Cependant dans un troisième exemple, re¬ 
latif â Sardine Ha Jussieu [= Sardine!la giàbosa (Bleeker, 1849)], J^d 2 augmente à 
chaque itération. La méthode n'est donc pas fiable. 


NOUVELLE MÉTHODE 

On part d'une valeur arbitraire de L choisie légèrement supérieure à 
la plus grande valeur moyenne L t observée. On en déduit une première valeur de 
k y k° pente de la droite de régression (3). En posant K = log(l — L f /L w ) et X = i 
réquation (3) s'écrit Y — — k°X -h k°t 0 * 

La_pcnte k° correspondant à est égale â eov(X t Y)jvarX. 

Or on a X — (/x + I )/ 2 varX = (rt 2 — 1 )/12 

et V = \ £log(l - U LL) 

On en déduit : 

=7^7 [Si fcwo - vo - ^r~ - vo] 

qui peut s'écrire: 


k° = cY togu - i,/o + dYt io g (i - 4/o (5) 

en posantc = _i_ etrf = __ii_ 

Cette valeur de k° est fonction de ZO seulement, k° — g(L^). 

Si l'on porte celte valeur de k* dans la relation (4), on obtient une nouvelle 
valeur de Ll 3 : 

k Q 

l l _ "V - S 2 

(«-!)(« - O 

est donc fonction de soit H cette fonction, on peut écrire: 

LL - H{k°) = HUiLl)-] 


D'une façon générale: 


L 


r+ 1 

oc 


(n - !)(«*'- 1} 


( 6 ) 


avec k r = g{LQ |d après (5)]. 

On voit donc que 1 — L'expression analytique de ^{L M ) s'obtient 

immédiatement en remplaçant dans (6} k r par son expression (5), 

Si r tend vers l infmi, et L^ x tendent vers et le problème est de 
résoudre l'équation x — \î/(x) que l'on écrira f(x) — ij/(x) — x = 0, 
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La résolution de cette équation est obtenue par la méthode classique 
d'approximations successives par tes tangentes* On part de la valeur arbitraire 
— L» et on écrit réquadon de la tangenLe en ce point. 


- —f(x 0 ) l'intersection de cette droite avec l'axe des x a pour 

A*o) 


coordonnées (y = 0, Xj) avec x Y - x 0 — 


/(*.) 


En répétant le processus on a d'une façon générale: 

f( X r) 


x r+l =x„- 


f (*,) 


et dans le cas qui nous intéresse où /{x ) — — x: 


LL +l 



*(Q ~ C 

mrj - 1 


( 7 ) 


Avec cette méthode, il est inutile de calculer les valeurs intermédiaires du 
paramétre k à chaque itération* Cest seulement lorsque la valeur de est 
stabilisée, au degré de précision souhaité, que l ont peut facilement calculer les va¬ 
leurs correspondantes de k s t D et Yjd 2 * 

Nous avons testé la valeur de cette méthode en prenant les données de Ra- 
fail (1973) relatives â Sardine lia Jussieu* puis nous Lavons appliquée aux données 
originales que Lun de nous avait récoltées, relatives aux Trachurus trachurus de 
Tunisie* 


Première application : Sardineüa jussicu 

La valeur de a été choisie égale à 19 cm, la plus grande valeur observée 
pour L étant IS cm. 


Itération 1 


Itération 2 


Age 

L, 

I - 1^ L oo 

Log(l - Lj Lso) 

1 

n,oo 

0*421 

-0,S65 

2 

14,00 

0,263 

-1,335 

3 

15,50 

0,184 

-1,692 

4 

17,00 

0,105 

-2,151 

5 

18,00 

0,053 

-2*944 

00; k = 0,508; l = - 

0,581; X d 3 - 0,488. 

Age 

L, 

1 - L, Loo 

Lûg( I - L^Loo) 

1 

11,00 

0,440 

'0,822 

2 

14,00 

0,237 

-1,249 

3 

15,50 

0,210 

-1,559 

4 

17,00 

0,134 

-2,108 

5 

18,00 

0,033 

-2,490 


L m « 19,63; k = 0,410; t G - -0,969; £ d 2 = 0,134. 
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Itération 3 


Itération 5 


Age 

L, 

1 - L, Loo 

Log( 1 - L, Loo) 

î 

11,00 

0,445 

-0,809 

2 

14,00 

0,294 

-1,224 

3 

15,50 

0,21 g 

1,522 

4 

17,00 

0,143 

-1,947 

5 

18,00 

0,092 

-2,383 

83; k = 0,387 ; t„ = 

-1,074; X d 2 = 0.105. 

Age 

L, 

t-4 ■ u 

Log( 1 - 1^ Loo) 

1 

11,00 

0,446 

-0,808 

2 

14,00 

0,294 

-1,223 

3 

1 5,50 

0,219 

-1,520 

4 

17,00 

0,143 

-1,944 

S 

18,00 

0,093 

-2,378 

84; k = 0,386; l 0 = 

-1.080; Xd J = 0,103. 

Age 

Lt 

1 - L, / Loo 

Log( 1 - Loo) 

1 

11,00 

0,446 

-0,808 

2 

14,00 

0,294 

-1,223 

3 

15,50 

0,219 

-1,520 

4 

17,00 

0,143 

-1,944 

5 

18,00 

0,093 

-2,378 


« 19,84; k - 0,386; t 0 = -1,080; Vd 2 - 0,103, 


Puisque X d z ne diminue plus, il est inutile de poursuivre les itérations* 
Celles-ci auraient pu être arrêtées à la quatrième. 1 e modèle auquel on aboutit fi* 
nalement donne les estimations suivantes aux différents âges: 


Age 

(année) 

L observée 
(cm) 

L calculée 
(cm) 

î 

11,00 

10,94 

2 

14,00 

13,80 

3 

15,50 

15,73 

4 

17,00 

17,05 

5 

18,00 

17,94 


Pi en que Rafail se soit arrêté à la deuxième itération sa méthode donne à ta 
quatrième itération des valeurs proches de celles que l'on a trouvées, mais la 
somme des carres des écarts V d 2 augmente progressivement au fur et à mesure 
qu'on avance dans les itérations; elle passe de 0,153 â la première itération à 0,241 
à la quatrième itération pour an oindre 1,313 à la %ingtième itération, etc, 11 n'y a 
donc pas do convergence de JT d 2 et on ne sait pas à quel niveau il faut arrêter le 
calcul* 
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Deuxième application : Trachurus trachurus 

La valeur iniliale de a été choisie égale à 36,62 cm, la plus grande val¬ 
eur observée pour L t ayant été de 35,62 cm. 

Itération 1 


Age 

u 

1 - 1^ , Lao 

Log( 1 - 1^ Loo) 

1 

11,25 

0,693 

-0,367 

2 

16,17 

0,558 

-0,583 

3 

21,10 

0,424 

-0,858 

4 

25,03 

0,316 

-1,150 

S 

28,37 

0,225 

-1,490 

6 

30,81 

0,159 

-1,841 

7 

33,00 

0,099 

-2,314 

8 

34,56 

0,056 

-2,878 

9 

35,62 

0,027 

-3,601 


l-oo « 36,62; k = 0,390; t 0 = Ü/707; £ d 2 ^ 62J 24. 

Itération 2 

= 39,19; k = 0,264; t 0 = -0,222; V d 2 = 4,740 
Itération 3 

= 42,31; k = 0,197; t Q = -0,542; £ d" = 0,529 
ï te ration 4 

L w = 42,30; k = 0,197; t 0 = -0,542; £ â 2 = 0,529. 

Calcul arrêté à la 4ème itération. 

Puisque d 2 ne diminue plus, il est inutile de poursuivre les itérations, 
celles-ci auraient pu être arrêtées à la troisième. Le modelé auquel on aboutit fi¬ 
nalement donne les estimations suivantes aux differents âges; 


Age 

(année) 

L observée 
(cm) 

L calculée 
(cm) 

l 

11,25 

11,10 

2 

16,17 

16,69 

3 

21,10 

21,28 

4 

25,03 

25,05 

5 

28,37 

28,14 

6 

30,81 

30,68 

7 

33,00 

32,76 

S 

34,56 

34,47 

9 

35,62 

35,88 


CONCLUSION 

La méthode qui vient d'ètre exposée est simple et rapide car il suffit en 
général de trois à quatre itérations successives. Llle donne de meilleurs résultats 
que la méthode de Rafail e! surtout, elle est beaucoup plus fiable car elle 
s'applique à tous les cas. 

Les calculs ont été programmes au centre de calcul El Khawarezmï de la 
Faculté des Sciences de Tunis, où je programme est disponible (bande SG0237). 
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